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Chapter 2 


Complex Representation of 
Continuous-Time Periodic Signals 


Leonhard Euler 
1707 - 1783 


Leonhard Euler was a pioneering Swiss mathematician and physicist. He made important 
discoveries in fields as diverse as infinitesimal calculus and graph theory. He also introduced 
much of the modern mathematical terminology and notation, particularly for mathematical 
analysis, such as the notion of a mathematical function. Euler is considered to be one 
of the greatest mathematicians to have ever lived. A student of Johann Bernoulli, Euler 
was the foremost scientist of his day. Born in Switzerland, he spent his later years at the 
University of St. Petersburg in Russia. He perfected plane and solid geometry, created the 
first comprehensive approach to complex numbers. Euler was the first to introduce the 
concept of log x and e* as functions and it was his efforts that made the use of e,1 and 7 
the common language of mathematics. Among his other contributions were the consistent 
use of the trigonometric sine, and cosine functions and the use of a symbol for summation. 
A father of 13 children, he was a prolific man in all aspects, languages, medicine, botany, 
geography and physical sciences and has left his mark on our scientific thinking.— From 
Wikipedia 
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Euler’s Equation 


eft — coswt + jsinwt 


This equation is called the Euler’s equation. Bertrand Russell and Richard Feynman both 
gave this equation plentiful praise with words such as “the most beautiful, profound and 
subtle expression in mathematics” and “the most amazing equation in all of mathematics.” 
This perplexing equation was developed by Euler (pronounced Oiler) in the early 1800s. 


The e/®« in Euler’s equation is a decidedly confusing concept. What exactly is the role 
of j in e/®'? We know from algebra that it stands for /—1 but what is it doing here with 
the sine and cosine? Can we even visualize this function? 


The complex exponential 


The function e/° goes by the name of complex exponential (CE). This function is of the 
greatest importance in signal processing and Fourier analysis. We are going to discuss 
its conceptual nature and its application to Fourier analysis. 


et — coswt + jsinwt (2.1) 


Looking at Eq. (2.1), we see the complex exponential on the left side. It is called the 
positive complex exponential, for the simple reason that the exponent of e is positive. 
On the right of the equal sign, underlined, is the expression containing a sine and a 
cosine. For now, ignore the complex exponential e/“‘ on the left-hand side and examine 
only the right-hand side of this equation, containing the sine and cosine waves, with the 
complex operator j thrown in. 


This is a complex function which often means that it is a 3D function. It can be plotted 
by assigning a particular value to the radial frequency w, and then for a range of time t, 
calculating both cos wt and sin wt values. Because the value of w is held constant, we 
have three values, the independent time variable, t, and associated sin and cos values, 
cos(wt) and sin(wt). With these three values, t, cos(wt) and sin(wt), we can create 
the 3D plot shown in Fig. 2.1. Time is plotted on the x-axis, and y and z axis contain 
the values of the function cos(wt) and sin(wt). Each plotted helix is for a particular 
value of frequency, w. Changing the frequency will change this figure. 


The expression for the negative complex exponential (negative exponent) is written 
similar to Eq. (2.1) as: 


e Jt — cos wt —jsinwt (2.2) 


We can plot these figures in Matlab using this code. 
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(a) elt - (b) e it 


sin(wt) 
—sin(wt) 


Figure 2.1: e/®' when plotted looks like is a helix. It is a 3-D function of three values, time t, 
sin(wt) and cos(wt) for a fixed frequency w. The exponent of the CE indicates direction of 
advance or movement, (a) positive exponent and (b) negative exponent. 


(a) e" =cos(w,t) + jsin(w,t) (b) e 7 =cos(u,t) — jsin(w,t) 


Sine phase - 
is opposite. "| 


sin(w,t) 
-sin(w,t) 


Figure 2.2: The projections of the complex exponential are sinusoids. (a) e)°' and its two 
projections, (b) e/°* and its two projections. Note that the sine wave in (b) has different phase 
than one for the positive CE in (a). That is the only difference between the two CEs. 


t = 0:0.01:5; 

y=exp ((j*(2*pi))*t); %& positive exponent 
subplot (1,2,1); 

plot3(t,real(y),imag(y), 'linewidth', 1); 
subplot (1,2,2); 

y2=exp ((-j* (2*pi))*«t); & negativ xponent 
plot3(t,real(y2),imag(y2), '‘linewidth', 1); 


N DOD Ww fF Ww NO 


Projections of a complex exponential 


As the CE is a complex function, we examine its projections on the real and the imaginary 
axes. In Fig. 2.2(a) we plot the projections of the helix on the Real and the Imagi- 
nary planes. In Cartesian terms, these would be called the (X,Y) and (Z,Y) planes. 
The projections of the complex exponential on these two planes are sine and cosine 
waves. The Real projection of the complex exponential is a cosine wave and Imaginary 
projection is a sine wave. 


41 


CHAPTER 2. COMPLEX REPRESENTATION OF CT PERIODIC SIGNALS 


These projections can be plotted using the following Matlab code. 


ie = 020.0125> 
2 y=exp((j*(2*pi))«*«t); & Positive exponential 
3 subplot (1,2,1 
lot3(t,real(y),imag(y), 'linewidth', 1); 
h2 = linspace(-2,-2, length(real(y))); hold on; 
(t, real(y), th2,'linewidth', .5); hold on; 
h = linspace(2,2, length(real(y))) 
(t, th, imag(y), 'linewidth', .5) 
E(1,2;,2)3 
wo y2=exp((-—j*(2*pi))*t); © Negativ xponential 
1 plot3(t,real(y2),imag(y2), '‘linewidth', 1); 


~ 


wz th2 = linspace(-2,-2, length(real(y2))); hold on; 
133 plot3(t, real(y2), th2,'linewidth', .5); hold on; 
w4 th = linspace(2,2, length(real(y2))) 

i plot3(t,th, imag(y2), 'linewidth', .5) 


For the negative exponent, or the so-called negative complex exponential, the projection, 
sine wave, is flipped 180° degrees from the positive exponential, as shown in Fig. 2.2(a). 
Often this exponential is referred to as having a negative frequency; however, it is not 
really the frequency that is negative. From the definition of the negative exponent 
exponential in Eq. (2.2), we see that negative sign of the exponential results in the 
imaginary projection, the sine wave doing a 180° phase change, or equivalently being 
multiplied by —1. 


The real part of the positive CE as well as the negative CE is a positive cosine wave. 
Re(e/*) = cos wt Re(e J") = cos wt 


The imaginary part of the positive CE is a positive sine but is negative for the negative 
CE. 


Im(e/*')=sinwt  Im(e/°!) =—sinwt 
While the negative exponential has as its imaginary part a negative sine wave, the 
positive exponential has a positive sine as its imaginary part. The real part, which is a 


cosine, is same for both. We do not see any negative frequencies here, an idea generally 
associated with the negative complex exponential. 


Adding and subtracting the complex exponentials, e/°‘ and e~/*, and then after a little 


rearrangement, we get these new ways of expressing a sine and a cosine. 


* (etiot +e Jer) = + (cos wt + jsinwt + cos wt — jsinwt) 
2 2j (2.3) 


= cos(wt) 
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(a) cos(wpt) = de +e) (b) sin(ugt) = z(e —e**") 


0 


Figure 2.3: (a) Plotting (e+)? + eJ®*)/2 gives a cosine wave with zero projection on the 
imaginary plane (b) Plotting (e*)°' — e°')/2j gives us a sine wave with zero projection on the 
real axis. 


1 ; ; 1 

—(eror —e Jer) = —(coswt + j sin wt —coswt + jsin wt) 
2j 2j (2.4) 
= sin(wt) 


Let’s see graphically what Eq. (2.3) and Eq. (2.4) look like. When we plot the two 
composite exponentials, we get the two plots as shown in Fig. 2.3. Fig. 2.3(a) shows 
that this composite exponential has a real projection of a cosine and Fig. 2.3(b), only 
the sine. The helix is gone, it has collapsed into a cosine or a sine. Hence, the sine and 
cosine can be said to be composed of complex exponentials. So a CE is composed of a 
cosine and a sine and conversely a sinusoid is also composed of two CE. 


We are so used to thinking of sine and cosine as sort of atomic functions, divinely 
given that it seems hard to believe that they can be created from other functions. But 
Eq. (2.3) and Eq. (2.4) tell us that both sine and cosine can be created by adding 
complex exponentials. This is a case of two 3D functions coming together to create 
a 2D sinusoid.This sounds strange but it’s actually not an unfamiliar concept. We can 
add two 2D functions and get a 1D function. An example is when we add a sine and a 
180° shifted sine, we get a straight line, a 1D function. So a 2D function created by two 
3D functions should not be a big stumbling block. 


The Sinusoids 


So how did Euler’s equation come about and why is it so important to signal processing? 
We will try to answer that by first looking at Taylor series representations of the expo- 
nential e*, sines, and the cosines. The Taylor series expansion for the two sinusoids is 
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Amplitude 


Figure 2.4: Sine wave as a sum of many exponentials of different weights. 


given in 2D by the infinite series as: 


x? x* x 
OS ay er 
: : : 2.5 
. x3 x x? : : 
snx=x—-—+—- 2+ 
3! 5! 7! 


Note that each one of these series is composed of many individual exponential functions. 
Thus, sine waves really are composed of exponentials! However, these are real expo- 
nentials that are nonperiodic and not the same thing as the complex exponentials. 
Complex exponentials are a special class of real exponentials and are used as an alternate 
to the sinusoids in Fourier analysis, as they are periodic and offer ease of expression and 
calculation, which is not obvious when we encounter them at first. 


Real exponentials are used in Laplace analysis as the basis set, instead of the complex 
exponentials we use in Fourier analysis. Real exponentials are more general than complex 
exponentials and allow analysis of nonperiodic and transient signals. Laplace analysis 
is a general case of which Fourier is a special case applicable only to periodic or mostly 
periodic signals. 


The Taylor series expansion for a real exponential e~ gives this expression 


te e., (2.6) 

2! 3! SO! 
Both of these equations, Eq. (2.5) and Eq. (2.6) are straightforward concepts after you 
have digested them a bit. And indeed, if these functions were plotted, we would get 
just what we are expecting, the exponential of e and the sinusoids. How close our plots 
come to the continuous function depends on the number of terms that are included in 


the summation. 


The similarity between the exponential and the sinusoids series in Eq. (2.5) and Eq. (2.6) 
shows clearly that there is a relationship here. Now let us change the exponent in 
Eq. (2.6) from x to j@. Note that the term @ is used here instead of wt simply to keep 
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the equation concise. Now we have by simple substitution, the e/° series: 


i iQ) iQ) iQ)4 95 
9 =14j04U PE VON MO Uys 
2! 3! 4! 5! 


(2.7) 


We know that j? =—1 and j#=1, 7 =—l, etc., substituting these values, we rewrite 
this series as 
: ; Q2 193 g4 iQ° g® iQ7 
Peng ee |: nla Sr ape e ALaeecal ee Le (2.8) 
2! 3! 4! S! 4! 7! 
We can separate out every other term with j as a coefficient to create a two-part series, 
one without the j and the other with 


62 694 96 
Ee ea a This is cosine 
3 395 +a7 (2.9) 
i oe" ae” ae er 
eae a This is j times sine 


It can be seen that the first part of the series is a cosine per Eq. (2.5) and the second part 
with j as its coefficient is the series for a sine wave. Hence, we see that Euler’s equation 
is really quite understandable! It came from the Taylor series. The concept of a Taylor 
series was formulated by the Scottish mathematician James Gregory but introduced into 
the world by the English mathematician Brook Taylor in 1715, so he gets his name on 
it instead. 


eJ® — cos(O) + j sin() 


We can now derive some interesting results from the Euler’s equation, such as the 
following. By setting 0 = 77/2, we can show that 


el"/2 = costa fZ) + j sin(/2) 
=0+j-1 
= 
By setting 9 = 37/2, we can show that 
gene = cos(3772) + j sin(37/2) 
=0+j--l 
=-j 


And an another interesting result: 


e/” = cos() + jsin(z) =—1 
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From this, we can write this amazing looking equation! 
e™+1=0 


So this complex looking function is not so complicated after all. These are some of the 
useful properties of a CE that we will be using. The function maintains a wondrous and 
mysterious quality, with added tinge of fear for some. However, we need to get over 
our fear of this equation and learn to love it. But, the question now is why bring up the 
Euler’s equation in context of Fourier analysis at all? Why all this rigmarole about the 
complex exponential, why are sines and cosines not good enough? They are certainly 
easier to visualize. 


In Fourier analysis, we computed the coefficients of sines and cosines (the harmonics) 
separately. In Chapter 1, we discussed the three different formulations of the Fourier 
series using sines and cosines, only with cosines, and with complex exponentials. Fourier 
analysis using the trigonometric form is not easy in practice. Trig functions are easy to 
understand but hard to manipulate. In fact, adding and multiplying them is a pain. On 
the other hand, doing math with exponentials is considerably easier. (See examples in 
Appendix A.) 


Using a single exponential representing both sinusoids can simplify calculations in Fourier 
series. This is the main advantage of switching to complex exponentials in using the 
complex form of the Fourier series. The math looks hard but is actually easier. However, 
complex exponentials bring with them some conceptual difficulties, which is that they 
are hard to visualize and are confusing at first. 


Typically, when we decompose something, we do it into a simpler form but here seem- 
ingly a more complex form is being employed. A simpler quantity, a cosine wave is now 
decomposed into two complex functions. However, the net result is that it will make 
Fourier analysis simpler. We will go from simplicity to complexity and then to simplicity 
again. 


Let us take this sinusoid that has a phase term to complicate things and present it in 
complex form. 


x(t) =Acos(wt + @) 
A 
=-e 
2 
— A jotyie 1 A,-jot.—j6 
2 2 


j(ot+6) 4 A ,-i(wt+6) 
2 
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In the last row, we separate the exponential into its powers. If we expand this expression 
into trigonometric domain using Euler’s equation, we see that indeed we do get back 
the trigonometric cosine wave we started with. 
A oh he 
a 5 (cosevt + 0) + jsinfeot +9) + cos(wt + 0) — jsineot+0)) 
=Acos(wt + 0) 


Fourier Series Representation using Complex Exponentials 


In Chapter 1, we used trigonometric harmonics (the sine and cosine) as a basis set to 
develop the Fourier series representation. The target signal was “mapped” onto a k set 
of sinusoidal harmonics, such as these based on fundamental frequency of wo. 


S= [ sin wot, COS Wot, SINZWot, cos2wot,...,coskwot,..., | 


A set of complex exponentials given by set S, can then be used alternately as a basis set 
for creating a complex (but preferred) form of the Fourier series. 


SS [ange Ot ee et Oe cae GOt | 
The big difference between these two forms is that whereas for the trigonometric form, 
the index k is strictly positive, for the complex exponential form, it covers both positive 
and negative integers. This is an important difference as we shall see in this chapter. 
These complex exponentials are an orthogonal set, making them easy to separate from 
each other. This is the main reason why we pick orthogonal signals to represent some- 
thing. Just as our 3D world is defined along three orthogonal axes, X, Y, and Z, 
our signals can be similarly projected on a K dimensional orthogonal set. We can not 
imagine these K orthogonal axes, but we know mathematically that they exist. 


As mentioned earlier, the Fourier series is a sum of weighted sinusoids. By weighted 
we mean that each sinusoid has its own amplitude and starting phase. The time is 
continuous but frequency resolution is not, frequency takes on discrete harmonic values. 
If the fundamental frequency is wa, then each w, is an integer multiple of wp and is 
discrete no matter how large k gets. The “distance” between each harmonic remains 
the same, wo. 


Let us examine the trigonometric Fourier series equation again. 
K K 
f(t)=ant > ay cos(wxt) + >) by sin(w,t) (2.10) 
k=1 k=1 
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The coefficients ag, a;, b, (which we call the trigonometric coefficients) are given by 
(from Chapter 1): 


Ti 


1 0 
dy = — f(t) dt 
0 Tete 
To 
a, = kaa f(t) cos(kwot) dt (2.11) 
To Jo 


aa tes 
b= = i f(t)sin(kwot) dt 
To Jo 


The presence of an integral tells us that time is continuous. Now substitute Eq. (2.3) 
and Eq. (2.4), as the definition of sine and cosine into Eq. (2.10), to get 


K K 
a . ; b ; . 
f(t) =ag+ » set + eee) + pS ae — e Tkwot) Q:12) 
k=1 k= “J 
Rearranging Eq. (2.12) so that each complex exponential is separated, we get 
1 1 
f(t)=ao+>) 5 (ae —JbneOt + yy, 5 (ak + ibe Teo" (2.13) 
k=1 k=1 


Let us redefine two new terms, called the complex coefficients, A; and By, as: 


1 . 

A, = 5 (tk =Jd_) (2.14) 
1 : 

By, = 3 (atk + jb,) (2.15) 


Substituting these new coefficients into Eq. (2.13), we get this representation. 
co co 
fo dg + > Apelkoot +) Bye Tkeot (2.16) 
k=1 k=1 


By substituting the complex exponential for the sinusoids, we can rewrite the trigonometric 
coefficients from Eq. (2.11) as: 


Ti 
2 = 1, ; 
nar i F(E\S (ert + eet) et ea) 
veer 2 


2” 1 
i t)—(elkeot —eJhoot) de 2.18 
k al Et Jozle e ) (2.18) 
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These can be expanded as follows. 


i are a i 
a=] f(telkorde+— | f(t)e To! de (2.19) 
k 
To 0 To (0) 
1 To 1 To 
be= S| fle deS =| f(cje Mort de (2.20) 
To 0 To (0) 


You can see that the trigonometric coefficient is split into two parts now, one for each 
of the exponentials. With some manipulation and using some simple complex math, we 
can show that the new complex coefficients are given by: 


T, 


1 : . 

Ap=— | f(t)e7*#o8 dt (2.21) 
To 0 
a tees 

Be=— | f(tel*@ot de (2.22) 
To 0 


Hence, A, can be thought of as the coefficient of the positive exponential and B, the 
coefficient of the negative exponential, analogous to the trigonometric coefficients be- 
longing to the cosine and the sine harmonics, a, and b,. The equivalence between the 
complex coefficients and the trigonometric coefficients takes the exact same form as the 
Euler equation. We will see in the example problems the effect of the trigonometric 
coefficients splitting into two parts. The total magnitude at a particular frequency is the 
same for both complex and the trigonometric forms, but appears split in the spectrum 
because of the sign of the index. We still have to contend with the term ay, the DC 
component in Eq. (2.16). We can get rid of this DC term by including it inside the 
summation. We do that by expanding the range of the index k from 0 to ©¢ instead of 
starting at 1. Rewrite Eq. (2.16) as 


[eoe) [o.e) 
F(t) =D iAgel*oo' +S Bye Tkeot (2.23) 
k=0 k=0 


Going forward, Eq. (2.23) can be simplified even further by extending the range of 
coefficients from —oo to co. Since the index now expands from positive domain to 
negative, we no longer need the negative CE in Eq. (2.23). Now both terms can be 
combined into one positive CE with a two-sided index so that we get a much more 
compact and elegant equation for the Fourier series. There is no negative exponential in 
this equation because the index takes care of that. And here is a much shorter equation 
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for the Fourier series in the complex domain. 


f= Ger (2.24) 


k=—oo 


The coefficient C; in Eq. (2.24) is given by 


F(tje Theo de (2.25) 


The equations Eq. (2.25) and Eq. (2.24) are called the complex form of the Fourier 
series. They are rigorously related to the trigonometric form. The magnitude calculated 
using the trigonometric form is exactly the same as the magnitude from this form. These 
two equations are the most used form of the Fourier series. The complex coefficients, 
C, are of course equivalent to the trigonometric coefficients by this relationship. 


1 
Ci =5 (a — jb;) ifk>0 
1 (2.26) 


As explained in Chapter 1, for the trigonometric form, the index k is always positive and 
therefore the spectrum for the Fourier series using the trigonometric form is one-sided. 
The x-axis for the one-sided spectrum is plotted against frequency starting at zero 
frequency or the “positive” fundamental frequency, which means that all k integer multi- 
ples of the fundamental frequency are positive. Because the index is positive, all frequen- 
cies are said to be positive. Now the index is allowed to be negative and this gives rise 
to the idea that the frequency has become negative, an effect of doing math in the 
complex domain. If the two forms are equivalent, does the frequency actually become 
negative when we use complex exponentials? With complex formulation, the index k 
spans from —©o to +00. We start with the negative index, go through calculations of all 
negative exponent exponentials and then the positive ones. Note that at no point is the 
fundamental frequency ever negative. Hence it is not the frequency of the exponential 
that is negative but just the index. The exponential with the negative index k is different 
from the positive exponential in that the sign of the imaginary part is negative. There 
is nothing here that implies that the frequency itself is negative. 


etIkeot — cos(kwot) + j sin(kwot) 


so (2.27) 


= cos(kwot) — j sin(kaw ot) 


We now equate this form with the trigonometric form that seemingly had only positive 
frequencies. To accomplish this, we look at what it takes to represent a sine and a cosine 
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using complex exponentials. 


1,; : 
cos(kwot) = 5 es e- skeoot) 
. Deets ae (2.28) 
sin(kwot) = ae ot _ pike ) 

J 


We require both a negative-index exponential and a positive-index exponential to des- 
cribe both the sines and cosines. Where index k is always positive on the left-hand 
side of this equation, it is both negative and positive on the right side. This traps us 
into thinking that frequency has changed sign. While in trigonometric form a positive 
index is enough to fully and completely represent the signal, in complex form it takes a 
double-sided index. The spectrum is plotting the product of the index and the frequency, 
(ka), and not just the frequency, w, on the x-axis. But we very quickly lose sight of this 
fact. We start talking about positive and negative frequencies because we confuse the 
range of the index with the sense of the frequency. 


Double-sided spectrum 


This gives rise to the double-sided spectrum which spans from values of k from -co to 
+oo. In this type of spectrum, the x-axis represents frequency, but that isn’t necessarily 
so. What we actually plot is the product of the harmonic index and the frequency. 
Calling it frequency gives us some intuitive comfort but then we have to worry about 
what a negative frequency means. When we say negative frequency, we have, in fact, 
unknowingly converted a complex idea into simple everyday language. Because of the 
plotting convention, the sign of the index is often forgotten and the axis is referred 
to as the frequency axis, spanning both positive and negative domains. In this book, 
we maintain that there is no such thing as a negative frequency. The idea comes from 
confusion caused by what the x-axis represents in the complex domain. 


The complex coefficient values are one-half of what they are calculated in the trigono- 

metric domain. Sometimes students think that this is because the frequency is being split 

into two parts, a negative part and positive part with each getting half the coefficient. 

This is how some authors try to explain the conundrum of positive and negative frequency 
in relation to Fourier analysis. However, they are just trying to explain a plotting con- 

vention. The real story is that we are not plotting frequency on the x-axis of a spectrum 

but the term +kw . First we call this "frequency" and then we try to explain this error. 

ACE is a helical function with a sense of rotational direction in addition to its frequency. 

The index k indicates this direction. 


The reason the amplitude values are split in half can be explained intuitively. We have 
let the index k go from —oo to +00, thus, each frequency is now multiplied by both 
a positive k and a negative k. However, in reality, each frequency has only a certain 
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amount of finite energy, so to make it all work out, the energy of this frequency is split 
between these two indices. Each index gets half. 


There are certain things that are defined only as positive quantities, such as volume, 
mass, age, etc., and frequency as a physical property is one of those things. It is never 
negative. In complex domain, we look at the signal in three dimensions, and hence a 
signal has a frequency and a rotational direction. So the term frequency, although a 
scalar, may be thought to include this other rotational parameter. 


Example 2.1. Compute the complex coefficients of a cosine wave. 


f(t) =Acos(wot) 


_ SZ eitkedent + Feith Veo! (2.29) 


This example is so simple that we can easily deduce the trigonometric coefficients just by 
looking at the expression. In fact, the equation is itself a perfect Fourier representation. 
The first thing we do is to write the equation in terms of its Euler form using Eq. (2.28). 
This form has the two CEs with the appropriate coefficients. The complex coefficients of 
the sinusoids are of magnitude A/2 located at k = 1 and k = —1. We plot the coefficients, 
C,, in Fig. 2.5 as the single-sided spectrum as well as the exponential coefficients, C;., as 
the double-sided form. The x-axis variable is kw. Since, wo. is a constant, we are really 
plotting, k, the index. Note, it is not the frequency that is negative but the harmonic 
index k. However, in a typical plot, the x-axis is labeled as frequency. In these plots, we 
have labeled it specifically as what it is, the term kw. 


(a) One-sided spectrum (b) Double-sided spectrum 


Real 


A 
2 


kay kay 


Figure 2.5: The spectrum of Acos wt. The trigonometric form gives us a one-sided spectrum with 
one component located at w. The complex form shows two components of half the total amplitude 
in the real plane. The vertical axis is amplitude for both forms. The total energy in both forms is 
the same, despite being split in exponential form in two parts. 


Example 2.2. Compute complex coefficients of a sine wave 


f(t) =Asin wot 
— A gitk=Noot _ A pik=—Deot (2.30) 
2j 2j 
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(a) One-sided spectrum (b) Double-sided spectrum 
A 
kay ka 
0 1@ . 


Figure 2.6: The spectrum of Asin wt. The trigonometric form gives us a one-sided spectrum with 
one component located the w. The complex form shows two components of half the total 
amplitude of opposite signs in the imaginary plane. 


This example is just the same as the cosine example. The single-sided spectrum is easy. 
It is simply a harmonic of magnitude A and located at k = 1 just as it is for the cosine 
wave in Ex. 2.1. In Eq. (2.30), we write the complex form based on Eq. (2.28), with the 
amplitudes of the two complex exponentials as A/2 and —A/2 located at k = +1. While 
the amplitudes were positive for cosine in Ex. 2.1, here they are of opposite signs. 


However, there is also the dreaded j in the denominator of Eq. (2.30). What to do with 
this? The presence of j tells us that the coefficients are on the imaginary axis, so they 
are to be plotted on the imaginary plane, right-angle to the plane on which a cosine lies. 
Drawn in 2D form it has no computational effect, only that the horizontal axis is called 
the imaginary axis. However, when there are both cosines and sine waves present in 
a signal, the coefficients of these two have to be combined not linearly but as a vector 
sum as seen in Fig. 2.7. This is because the harmonics are orthogonal to each other. 
When plotting the magnitude, it no longer falls in just the Real or the Imaginary planes 
but somewhere in between. In the next example, we combine both of these sinusoids 
in a single function. 


Re 


From sines 
plotted on 


Im Imag. plane 


Figure 2.7: Magnitude of the resultant vector if a signal contains both a sine and a cosine. 


Example 2.3. Compute the coefficients of f(t) = A(cos wt +sinwt). 


We can write this wave as 
A. A; A. A _. 
f(t) = =f + eI! + eh — Se IM (2.31) 
2 2 2j 2j 
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Once again, we convert this function to its complex form using Eq. (2.28). From this 
form, we can pick out the trigonometric coefficients easily. It is simply A for the cosine 
and A for sine with magnitude equal to square root of /2A located at w = 1. 


The coefficients of the complex exponentials, too can be obtained by looking at the 
coefficients of the two exponentials in Eq. (2.31). The e/°* exponential has two coeffi- 
cients, at 90° to each other, each of magnitude A/2. The vector sum of these is A//2. 
Same for the negative exponential, except the amplitude contribution from the sine is 
negative. However, the vector sum or the magnitude for both is the same and always 
positive. This is shown in Fig. 2.8(c) drawn in a more conventional style showing 
only the vector sum. Note that the total energy (sum of magnitudes) of the one-sided 
spectrum is exactly the same as that of complex. 


(a) (b) 


A C, 
Cc A k 
k 2 of 
2A A * 
1 : a 
—-lo : 2 
> lo ke 
k = Y 
0 lo i 3 A 
iad 2 
@ [ 
A A 
V2 V2 
ko 
-lo 0 lo 


Figure 2.8: Amplitude spectrum of Asinwt + Acos wt. 


Example 2.4. Compute coefficients of the complex signal f(t) =A(cos wt + jsinwt). 


This function is different from the one in Ex. 2.3 in that the sin is located in the imaginary 
plane. This is a complex signal. We can write this in CE form as: 


f(= See +See +See seve = Aeloe 


The coefficients from the sine and the cosine for the negative exponential cancel each 
other. On the positive side, the contribution from sine and cosine are coincident and 
add. Therefore, we see a single value at the positive index of k = 1 only. For this signal, 
both the single-sided and double-sided spectrum are identical. This is a surprising and 
perhaps a counter-intuitive result. 


Important observation: Only real signals have symmetrical spectrum about the origin. 
Complex signals do not. 
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Ic,| 


ka 
-lo 0 lo 


Figure 2.9: Double-sided spectrum of Acos wt + jAsin wt. 


Example 2.5. Compute the coefficients of this signal, f(t) =A. 


We can write this function as an exponential of zero frequency. 


f (t) =Acos(w = 0)t 

~A,i(o=0t , A,-i(o=0)t 
2 2 

=A. 


The trigonometric coefficient is A at w = 0. For the complex representation, we get two 
complex coefficients, both of amplitude A/2 but both at k = 0 so their sum is A, that is 
exactly the same as in the trigonometric representation. The function f(t), a constant, 
is a non-changing function of time and we classify it as a DC signal. The DC component, 
if any, always shows up at the origin. The single-sided and double-sided spectrum here 
are same as well. This signal has only one component at w = 0; hence, it has only one, 
coefficient dg, all the others are zero. Therefore, that is what we are seeing in Fig. 2.10. 
Just the dg coefficient plotted. 


Important observation: A component at zero frequency means that the signal is not zero- 


mean. 


C 


k 


ka 
0 


Figure 2.10: Double-sided spectrum of a constant signal of amplitude A. A constant signal is same 
as a DC value, hence, its spectrum always appears as an impulse at the origin. 
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Example 2.6. Compute the coefficients of x(t) = 2 cos?(wt). 
We express this function in complex form as: 


ejet as e jet i: 


x(t)= 2{ ; 


=1+ 1 yjewe 4 1 -j20t 
2 2 


C, 


|e 
Re 


ko 


—20 re) 20 
Figure 2.11: Double-sided amplitude spectrum of 2.cos?(wt). 


The coefficients become obvious when we put the trigonometric form into the complex 
exponential form. There is just one frequency, 2w and it has an amplitude of 5 at k of 
+2. The coefficient of zero frequency is 1. 


Observation: A squared signal by definition is always positive so the spectrum has a zero- 
frequency component at the origin. 


Example 2.7. Compute the coefficients of x(t) = 2cos(wt) cos(2wt). 


We can express this signal in complex form by making use of this trigonometric identity: 
2(cos(a) cos(b)) = cos(a — b) + cos(a + b). 


x(t) = cos(wt) + cos(3wt) 
= 1 ojo +4 1 jet a 1 J j3et +4 1 |-j30t 
2 2 2 2 
Of course, doing this in trigonometric form would have been just as easy; however that 


is not always true. We draw the spectrum as in Fig. 2.12. 


-40 -30 -2@ -1lo O lo 20 30 4a 
Figure 2.12: Double-sided amplitude spectrum of 2.cos(wt)cos(2wt). 


This signal has four coefficients, a, and a3; all the others are zero. So, that is what 


we are seeing in Fig. 2.12. 
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Example 2.8. Compute the complex coefficients of this real signal. 
f(t) =sin(4zt)+ 0.8 cos(8zt) + 0.3 sin(147t) 


We want all these frequencies to fall on a harmonic, so we select as the fundamental 
frequency, f = 1. The signal has six components, at k = +2,+4,-+7. We can write this 
equation in complex form as: 


1 »4- 1 oe 
th= = elan(k=2)t _ = pl2m(k=—2)t 
FO=5 5 


4 08 ojantk=aye 0:8 jan(k=-4)t 


2j 2j 
4 03 ojamtk=7)t_ 9:3 jan(k=-7)t 
2j 2 


The contributions at k = 2 comes only from a sine, at k = 4 from a cosine and at k = 7 
only from a sine. Note, we plot these on the same line at full magnitude as if j does not 
exist in the equation. (We will drop mentioning the index k and call it frequency to fall 
in line with common usage. However, note that it is this sloppiness in terms that causes 
us to question our sanity and start asking: what is a negative frequency?) 


Ic,| 


4 4 


et 


-7o -4@ -20 0 20 40 70 


Figure 2.13: Two-sided symmetrical magnitude spectrum. 


This signal has six coefficients, a2, a4, and a,7. All the others are zero. So that is what 


we are seeing here in Fig. 2.13. Just those three coefficients plotted on each side, with 
half the amplitudes from the time domain equation. 


Example 2.9. Compute the complex coefficients of this real signal with phase terms. 
Then compute its power spectrum. 
x(t)=3+6cos(4nt + 2)+4jsin(4nt +3)—3jsin(107t +1.5) 
We convert this to the complex form as 
x(t)=3+ (3e/4"* ei? + Se ite) + (2e/47* e/8 — Qe eo?) 

_ (Seft0rt ets _ Se H10t g-j1.5) 

= 3+ el4"*(3e!2 + 2el9) + e 14" (Se I? — 2eI3) + e/1*(1 5e/15) 

rf eH10nt (7 Sei15) 
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The magnitudes of the exponentials come from the terms in parenthesis. To add, we 
need to convert them first to rectangular form as follows. (see Appendix A). The CE 
e/47 has the following coefficients: 


el4tt _, (3e!? + 2e/%) 


=> |[3e/2 4 2e/3| = (3 cos(2)+ 2 cos(3))” + (3 sin(2) +2 sin(3))” 
= 4.414 


Similarly, the coefficient of the negative exponential is 


e J4nt _, (3eJ? — 2e/) 


=> |3e 2 — 227) = V3 cos(2) —2 cos(3))° + (3 sin(2) + 2 sin(3))” 
= 2.552 
We draw the spectrum in Fig. 2.14 and note that the spectrum is not symmetric because 
the signal is complex. 


Important Observation: Most signals we work with are complex and their spectrum are 
rarely symmetrical. 


(a) Magnitude (b) Power, dB 
ic, 


1.5 


5a 20 0 20 “5a ko -5o0 --20 O 20 5a 


Figure 2.14: (a) Two-sided Magnitude spectrum of a complex signal, (b) its power spectrum 
computed by squaring each component and converting to dBs. 


Power spectrum 


You are familiar with this expression of power from circuits. The power is defined as: 
P=V2/R 


Here, V is the voltage or the amplitude of the signal and R the resistance. Assume that R 
is equal to 1.0; in this case, the normalized power is equal to the voltage squared. If we 
square the peak voltage, we get peak power, and if we square mean voltage, then mean 
power is obtained. This idea is exactly the same as Parseval’s theorem, which states 
that the power in a particular harmonic is equal to the square of its amplitude. So, for 
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this particular example, to obtain the power spectrum, we just square each magnitude, 
convert it into dBs, and then normalize for maximum power. The result is shown in 
Fig. 2.14(b). 


Now a difficult but an important example, a periodic signal of square pulses. 
Example 2.10. Compute the Fourier coefficients of a periodic signal consisting of square 


pulses. The square pulse is of amplitude 1v that lasts tT seconds and repeats every T 
seconds. (Note that t and T are different and independent quantities.) 


—1/2 7/2 
Figure 2.15: A square wave of period T and duty cycle 2/T. 


First, we compute the coefficients for a general case for pulse time equal to t seconds 
and repeat time, or the period of the wave, equal to T seconds. The term T/T is called 
the duty cycle of the wave. 


1 T/2 
CG. = al 1 e J2knfot dt 
—7/2 


Note that outside of —t/2 < t < t/2, the function is zero. This integral is given by 


1 eo Ui2nk/TI)t t/2 
eS a ap 
T —j2nk/T —</2 
1 (ei2tk/T)E — o(-J2nk/T)$ 1 
7 r( j2 ean 
_ tsin(knt/T) 
~ T kn(t/T) 


Replacing the duty cycle term tT with term r, the equation becomes easier to understand. 


in(k 
C, = ae =rsinc(kzr). 
kar 
Let us set the duty cycle of this signal to 0.5. 
T 1 
r=—=-= 
T 2 
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Amplitude 


Amplitude 


Amplitude 


Time, t 


Figure 2.16: Fourier coefficients of a square pulse. 
Note that as the pulse gets narrower, its frequency response gets shallower. However, as the duty 
cycle increases, such as in (c), the response is becoming narrower. 


Substituting r in the above equation, we can write the expression for the coefficients of 
this signal as: 


C. = ssine(kr/2) (2.32) 


This is the sinc function. It comes up so often in signal processing that it is probably the 
second most important equation in DSP after the Euler’s equation. 


Now, we can plot the coefficients of the repeating continuous-time square pulse coeffi- 
cients for various duty cycles. Note how the peakedness of the main lobe changes 
inversely with the duty cycle. A narrow pulse relative to the period in Fig. 2.16(a), 
has a shallower frequency response than one that takes up more of the period. The zero 
crossings occur at inverse of the duty cycle. For r = 0.5, the zero crossing occurs at 
k = 2, for r = 0.25, the crossing is at k = 4, and for r = 0.75, the crossing occurs at 
n= 1.33. At r = 1.0, the pulse would be a flat line and it will have an impulse as its 
frequency response. For very small r, the pulse is delta function-like and the response 
will go to a flat line. Note the usage of words, frequency response. This is just another 
name for the spectrum. 


Although the equation for this function is fairly easy, it takes a little while to develop 
intuitive feeling. We cannot overemphasize the importance of this signal and you ought 
to spend some time playing around with the parameters to understand the effect. We 
will of course keep coming back to it in the next few chapters. 


In this chapter, we covered the complex form of the Fourier series as a prelude to the 
next topic, the Fourier transform. We see that even though the time domain function is 
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continuous and periodic, the Fourier series coefficients (FSC) and hence the spectrum 


developed is discrete. In Chapter 3, we will see how that affects the Fourier analysis. 


Summary of Chapter 2 


In this chapter, we looked at the complex exponential as a concise way of representing 


the Fourier series equation. The complex exponentials make the Fourier series math 


easier. The spectrum however is now shown as double-sided, which means that the 
frequency index spans from —0co to +00. This has the effect of splitting the trigonometric 


coefficients into half. 


Terms used in this chapter: 


Euler’s equation 


Continuous-time complex exponential, e/°* 


of frequency w. 
Complex coefficients of the Fourier series, C;, 


Double-sided spectrum 


. The Euler equation defines a complex exponential as a 3D function consisting of 


a cosine and sine in quadrature. 


e/°' =coswt + jsinewt 


e /°' =coswt —jsinwt 


. The Euler expression of a CE can have a positive or a negative exponent. The sign 


change indicates a change of direction of the function. 


. The sine and cosine alternately can be represented by two CE of different signs. 


We use the following expressions to represent them in the Fourier series to obtain 
a complex form of the Fourier series equation. 


ais. 
cos(kwt) = ie + gue) 


jee 
sin(kwt) = 7 ae = esha) 


. To represent a periodic signal using the complex exponentials requires a double- 


sided harmonic index k, unlike the trigonometric case where the harmonic index 
is positive. 


. The harmonic index extends from —K <k <+K. K can span from -0o to +00. 
. The x-axis now represents values from —Kw, < kwo < +Kwo and this is often 


read as representing negative frequency when, in fact, it is the index that is 
negative. 
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7. The FSC instead of being of three types can now be represented by a single 


equation: 
To/2 
f (t)e Jk! de 
—T)/2 


1 


== 
k Ty 


8. The fundamental properties remain the same, the time in this representation is 
continuous and frequency is discrete with index k, which is an integer. 

9. The amplitude spectrum obtained from the complex representation looks different 
from the one-sided spectrum. We call this spectrum a two-sided spectrum. The 
amplitude value for a particular harmonic is now split in half for the positive index 
(so-called positive frequency) and half for the negative index (so-called negative 
frequency). The Oth component, however, remains the same. This conserves 
energy and makes both forms equivalent. 
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Questions 


1. 
. What is the expression for a sine in the complex form? 


What is a complex exponential? 


. If we look along this axis, the CE appears as a circle. What axis is that? Which 


parameter is not visible in this view? 


3 a : = «3 
. What is the value of e~/”, e/?7, and eJ4”™, 
. Why is this equation true? e/7+1=0 
. What is the difference between two complex exponential of the same exponent 


but different signs, such as e/°«' and e«*. If we add these two signals, what do 
we get in the complex domain? 


. What dimensional space is required to plot a complex exponential? 
. The term phasor is often used in relation to complex exponentials, what is it? 
. If you plot a sinusoid plus its shifted version, sin(2wt) + sin(2wt + ¢o), what is 


the phase of the new signal? 


. Represent this sinusoid as a complex exponential: cos (cot + F). 
. What is the relationship of the Taylor series to a sinusoid? 
. What is the advantage of using complex exponential as a basis set instead of sines 


and cosines? 


. Given these CEs, give their expression in the Euler form: e 3, (e 34 + e J), and 


ert, 


. How would you express phase in the complex exponential form? 
. Write the CE form of these signals: 


(a) sin(7zt + a) 
(b) cos(7xt + F)—sin(7nt — F) 
(c) cos(3nt + 5) 


. What is the magnitude of this complex signal, e~°/‘? 
. The real part of a CE has a peak amplitude of 1 and the imaginary part has the 


peak amplitude of 2. What is the peak and the average power of this signal? 


. How do we transmit a complex signal? What does it look like? 
. What does division by j mean? 

. What does multiplication by j mean? 

. What is the magnitude and phase of these signals: 


(a) ( sin(wt)— cos(ct)) 
(b) (5 cos(wt) — sin(wt)) 
(c) (2 cos(wt) — 5 sin(cwt)) 


. What is a single-sided spectrum? What does it represent? 
23: 
24. 


Given the amplitude spectrum, how would you compute the power spectrum? 
What is two-sided spectrum of these signals? 

(a) f(t) =—cos(2zt) + cos(9nt) + sin(127t) 

(b) f(t) =2cos(9rt) + 2cos(18z2t) 
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25. 
26. 


27. 


28. 
29. 
30. 
31. 
32. 


33. 


When plotting a two-sided spectrum, what does the x axis represent? 

If you are given the real and imaginary components of a signal, how do you 
compute the phase? Is phase changing with time or frequency? 

For a complex signal, both real and imaginary signals can have nonzero phase, so 
what is the phase of a complex signal? How is it different from the phases of the 
components? 

What is the relationship of the trigonometric coefficients to the complex coefficients? 
How do explain the idea of negative frequency? 

Why the spectrum of a complex signal is always one-sided? 

What is the mean and the peak power of a sinusoid of amplitude 1.0? 

If we have a periodic signal of square pulses with a duty cycle of 0.1. How much 
wider is its spectrum as compared to a pulse that has a duty cycle of 0.5? 

What happens to the magnitude spectrum if phase of the signal changes? 
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Appendix A: A little bit about complex numbers 


We can use complex numbers to denote quantities that have more than one parameter 
associated with them. A point in a plane is one example. It has a y coordinate and a 
x coordinate. Another example is a sine wave: it has a frequency and a phase. The 
two parts of a complex number are denoted by the terms Real and Imaginary, but the 
Imaginary part is just as real as the Real part. Both are equally important because they 
are needed to nail down a physical signal. 


The signals traveling through air are real signals and it is only the processing that is done 
in the complex domain. There is a very real analogy that will make this clear. When 
you hear a sound, the processing is done by our brains with two orthogonally placed 
receivers, the ears. The ears hear the sound with slightly different phase and time delay. 
The received signal by the two ears is different and from this our brains can derive fair 
amount of information about the direction, amplitude and frequency of the sound. So 
although yes, most signals are real, the processing is often done in complex plane if we 
are to drive maximum information. 


Figure A.1: Multiplication by j shifts the location Of a point on a plane by 90°. 


The concept of complex numbers starts with real numbers as a point on a line. Mullti- 
plication of a number by —1 rotates that point 180° about the origin on the number line. 
If a point is 3, then multiplication by —1 makes it —3 and it is now located 180° from 
+3 on the number line. Multiplication by —1 can be seen as 180° shift. Multiplying this 
rotated number again by —1, gives the original number back, which is to say by adding 
another 180° shift. Therefore, multiplication by (—1)* results in a 360° shift. What do 
we have to do to shift a number off the line, say by 90°? This is where j comes in. 
Multiply 3 by j, so it becomes 37. Where do we plot it now? Herein lies our answer to 
what multiplication with j does. Multiplication by j moves the point off the line. 


Question: What does division by j mean? 
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Figure A.2: a. A point is space on a Cartesian diagram. b. Plotting a complex function on a 
complex plane. 


Answer: It is same as multiplying by —j. 
ms x J = — = =Jx. 
J J 


This is essentially the concept of complex numbers. Complex numbers often perceived 
as “complicated numbers” follow all the common rules of mathematics. Perhaps, a 
better name for complex numbers would have been 2D numbers. To further complicate 
matters, the axes, which were called X and Y in Cartesian mathematics are now called 
Real and Imaginary, respectively. Why so? Is the quantity j3 any less real than 3? 
This semantic confusion is the unfortunate result of the naming convention of complex 
numbers and helps to make them confusing, complicated and, of course, complex. 


Now let us compare how a number is represented in the complex plane. Plot a complex 
number, 3+ j4. In a Cartesian plot, we have the usual X — Y axes and we write 
this number as (3,4) indicating 3 units on the X-axis and 4 units in the Y-axis. We 
can represent this number in a complex plane in two ways. One form is called the 
rectangular form and is given as: 

Z=xt+jy 


The part with the j is called the imaginary part (although of course it is a real number) 
and the one without is called the real part. Here 3 is the Real part of z and 4 is the 
Imaginary part. Both are real numbers of course. Note that when we refer to the 
imaginary part, we do not include j. The symbol j is there to remind you that this 
part (the imaginary part) lies on a different axis. Hence writing these as: Re(z) = x and 
Im(z) = y Alternate form of a complex number is the polar form: z = MZ0@ 


where M is the magnitude and @ its angle with the real axis. The polar form that looks 
like a vector and in essence it is, is called a Phasor in signal processing. This idea comes 
from circuit analysis. We also use it in signal processing but it seems to cause some 
conceptual difficulty. Mainly because, unlike in circuit analysis, in signal processing 
time is important. We are interested in signals in time domain and the phasor, which is 
a time-independent concept, is confusing. 
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Question: If z = Ae/®! then what is its rectangular form? 
Answer: z =Acos wt + jAsinwt. 


We just substituted the Euler’s equation for the complex exponential e/®'. Think of e/°¢ 
as a shorthand functional notation for the expression cos wt + jsinwt. The real and 


imaginary parts of z are given by 


Re(z) =Acoswt and Im(z) =Asinwt. 


Converting forms 
1. Given a rectangular form z = x + jy then its polar form is equal to 


ere J/x2+ y2Ztan! x ifx >0 
7 Vf x2 + y2Z(tan-1~ + 1) ifx <O 


2. Given a polar form M Z@ then its rectangular form is given by 
x+jy=Mcos@+jM sind 
Example 2.11. Convert z = 520.927 to rectangular form 
Re(z) = 5 cos(0.927) =3 


Im(z) = 5 sin(0.927) = 4 
>z=3+)j4 


Example 2.12. Convert z = —1—j to polar form 
M = y(-1)? + (-1)2 = V2 
6 = arctan ea +7, sincex <O 
x 
—1 
= arctan —+7=-7 
—1 4 
3 
>z= V2LoR. 


Example 2.13. Convert z = 1+ to polar form. 


M = V(1)?+(1)? = v2 


6 = arctan om since x > 0 
x 


1 1 
=arctan- =—-7 
1 4 


=>z=V2/n/4. 
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Adding and Multiplying 


Add in rectangular form, multiply in polar. Its easier this way. 


1. Given z; =a+jb and z,=c+jd then z, +z, =(at+c)+j(b+d). 
2. Given 24> M,2Z0, and 22> MyZ6, then 21°29 >= M,M,Z(6; + 05). 


Example 2.14. Add z, = /220.785 and z = 520.927. 
Convert both to rectangular form. 
2, =1+jandz,=3+j4 
> 23 =(1+3)+j+4)=44+ 5. 
Example 2.15. Multiply 2] = 1+j and z.=3+ j4. 


First convert to polar form and then multiply. Although multiplying these two complex 
numbers in rectangular format looks easy, in general that is not the case. Polar form is 
better for multiplication and division. 


21 +B = V2Z0.785 x 520.927 = 5V2/1.71 
1 2 


Example 2.16. Divide z; = 1+ j and z. =3+4 j4. 


z,  72Z0.785 5 


ee 0.140 
Zz, 520.927 V2 


Conjugation 


The conjugate for a complex number z, is given by z* = x—jy. For a complex exponential 
e/°* is the complex conjugate of e/®*. In polar format, the complex conjugate is same 
phasor but rotating in the opposite direction. 


Rule: Ifz =MZ0@, then z* =MZ— 8. 


Useful properties of complex conjugates 
|z|? = zz* 


This relationship is used to compute the power of the signal. The magnitude of the 
signal can be computed by half the sum of the signal and its complex conjugate. Note, 
the imaginary part cancels out in this sum. 


1 
|z| = iG +2"), 
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